基于蒙特卡洛模拟与PMF模型的黄河流域沉积物重金属污染评价及源解析

沉积物是河流的重要组成部分,而沉积物中重金属的富集严重威胁着水环境安全.黄河流域沿程分布着众多工业城市,且流域水土流失量大,泥沙携带重金属进入河流导致沉积物重金属污染问题日趋严重,研究黄河流域沉积物中重金属污染状况对流域生态安全具有重要意义.收集2000~2020年发表的关于黄河流域沉积物中重金属（铅、镉、铬、砷、锌、铜、镍和汞）含量的数据,首先基于描述性统计及地统计法分析重金属的空间分布特征,进一步采用蒙特卡洛法进行地累积指数（I_geo）、潜在生态风险及毒性单位概率的评价,最后结合正定矩阵因子分解模型（PMF）与Pearson相关性分析确定污染源个数及贡献率.结果发现,黄河流域沉积物中ω（Pb）、ω（As）、ω（Zn）、ω（Ni）、ω（Cu）、ω（Hg）、ω（Cr）和ω（Cd）的均值分别为26.92、11.78、87.17、31.13、24.96、0.07、73.36和0.58 mg·kg^-1,分别超过黄河流域各省土壤背景值均值1.27、1.08、1.26、1.05、1.09、2.32、1.14和5.95倍,其中Cd超标倍数最大,应当引起重视;I_geo：Cd>Hg>Cr>Cu>Pb>Zn>As>Ni,Cd和Hg存在中度-严重污染;数据表明黄河流域上、中和下游沉积物重度生态风险占比分别为18.6%、15.7%和7.1%,呈递减趋势;黄河流域沉积物中重金属处于低毒性状态;溯源分析表明黄河流域沉积物重金属的4个来源分别是矿业源（42.2%）、自然活动（38.3%）、农业活动（11.6%）和电镀废水（7.9%）.研究结果可为黄河流域制定相关污染防控措施提供依据.     

中国典型流域有机磷酸酯的污染特征与风险评估

有机磷酸酯（OPEs）作为溴代阻燃剂的替代品,在生产生活中被广泛使用,其环境污染和毒性效应受到广泛关注.梳理中国七大典型流域水体中OPEs的含量水平和分布特征,计算中国成人、青少年和儿童的OPEs日均饮水暴露量,对其健康风险进行评估,并通过蒙特卡洛模拟方法评价了结果的可靠性.检索整理了12种OPEs对水生生物的毒性效应浓度,构建物种敏感度分布（SSD）曲线,对生态风险进行评估.结果表明,低等暴露水平下,七大流域ΣOPEs的第5百分位浓度为52.61 ng·L^-1;中等暴露水平下,七大流域ΣOPEs的中位浓度为499.74 ng·L^-1,磷酸三氯乙酯（TCEP）、磷酸三乙酯（TEP）和磷酸三（1,3-二氯-2-丙基）酯（TDCP）为主要污染物;高等暴露水平下,七大流域ΣOPEs的第95百分位浓度为1904.4 ng·L^-1,是其中等暴露水平的3.8倍,长江流域的ΣOPEs浓度最高.健康风险评估表明,不同人群通过饮用水暴露于OPEs的非致癌风险均在可接受范围内.磷酸三甲酯（TMP）、磷酸三异丁酯（TiBP）和TCEP是致癌风险的主要贡献者.生态风险评估结果表明,TCEP在高等暴露水平下存在中等生态风险;磷酸三正丁酯（TnBP）在中等暴露水平下存在中等生态风险,在高等暴露水平下,有较高生态风险;磷酸三苯酯（TPhP）在低中高3种暴露水平下的风险商均大于1,有较高生态风险,需要重点关注.     

VOC emission caps constrained by air quality targets based on response surface model: A case study in the Pearl River Delta Region, China

Because of the recent growth in ground-level ozone and increased emission of volatile organic compounds (VOCs), VOC emission control has become a major concern in China. In response, emission caps to control VOC have been stipulated in recent policies, but few of them were constrained by the co-control target of PM_2.5 and ozone, and discussed the factor that influence the emission cap formulation. Herein, we proposed a framework for quantification of VOC emission caps constrained by targets for PM_2.5 and ozone via a new response surface modeling (RSM) technique, achieving 50% computational cost savings of the quantification. In the Pearl River Delta (PRD) region, the VOC emission caps constrained by air quality targets varied greatly with the NO_x emission reduction level. If control measures in the surrounding areas of the PRD region were not considered, there could be two feasible strategies for VOC emission caps to meet air quality targets (160 µg/m³ for the maximum 8-hr-average 90th-percentile (MDA8-90%) ozone and 25 µg/m³ for the annual average of PM_2.5): a moderate VOC emission cap with <20% NOx emission reductions or a notable VOC emission cap with >60% NO_x emission reductions. If the ozone concentration target were reduced to 155 µg/m³, deep NO_x emission reductions is the only feasible ozone control measure in PRD. Optimization of seasonal VOC emission caps based on the Monte Carlo simulation could allow us to gain higher ozone benefits or greater VOC emission reductions. If VOC emissions were further reduced in autumn, MDA8-90% ozone could be lowered by 0.3-1.5 µg/m³, equaling the ozone benefits of 10% VOC emission reduction measures. The method for VOC emission cap quantification and optimization proposed in this study could provide scientific guidance for coordinated control of regional PM_2.5 and O₃ pollution in China.     

Do microbial processes regulate the stability of a coral atoll's enclosed pelagic ecosystem?

Complex marine ecosystems contain multiple feedback cycles that can cause unexpected responses to perturbations. To better predict these responses, complicated models are increasingly being developed to enable the study of feedback cycles. However, the sparseness of ecological data often limits the direct empirical parameterization of all model parameters. Here we use a Bayesian inverse analysis approach to synthesize empirical data and ecological theory derived from published studies of a coral atoll's enclosed pelagic ecosystem (Takapoto Atoll, French Polynesia). We then use the estimates of flux magnitudes to parameterize probabilistic compartment models with two forms of heterotrophic consumption: (1) “bottom-up” donor-controlled heterotrophic consumption and (2) “top-down” mass-action heterotrophic consumption. We explore how the flux magnitudes affect the ecosystem's stability properties of resilience, reactivity, and resistance under both assumptions for heterotrophic consumption. The models suggest that the microbial uptake of dissolved organic carbon (DOC) regulates the long term rate of return to steady state following a temporary or pulse perturbation (resilience), and the cycling of carbon between abiotic pools and heterotrophic compartments regulates the short-term response (reactivity). In the bottom-up process model, the sensitivity of steady state masses following a sustained or press perturbation (resistance) is highest for the DOC pool following a sustained change to the microbial uptake rate of DOC. Further, a change in the microbial uptake of DOC propagates through the ecosystem and affects the steady state values of zooplankton. The analysis suggests that the food web is highly dependent on the recycling between the abiotic and biotic carbon pools, particularly as mediated by the microbial consumption of DOC, and this recycling determines how the ecosystem responds to perturbations.     

Proposing an information criterion for individual-based models developed in a pattern-oriented modelling framework

Individual-based models (IBMs) have been improved in quality and reliability in recent years with an approach called pattern-oriented modelling (POM). POM proposes guidelines to develop models reproducing multiple patterns observed on the field and to test systematically how well the IBMs reproduce them. POM studies used generally traditional methods of goodness of fit such as the sum of squares evaluation or ad hoc comparisons of fitting errors and variations. Model selection, however, can be a rigorous statistical approach based on information theory and information criteria such as the Akaike's information criterion (AIC) or the deviance information criterion (DIC). So far, it has not been tried to link POM to these rigorous techniques. The main problems to achieve that are: (a) the difficulty to have likelihood functions for IBMs’ parameters and (b) the possibility to obtain posterior distributions of IBMs’ parameters given the patterns to reproduce. In a first part, this paper answers problem (a) by proposing and explaining how to calculate a deviance measure (POMDEV) for models developed in a context of POM. And while answering the second problem, a second part of the paper proposes an information criterion for model selection in a POM context (the pattern-oriented modelling information criterion: POMIC). This criterion does not yet have the same theoretical foundation as, e.g., AIC, but uses formal analogies to the DIC. In a third part POMIC is tested with a modelling exercise. This exercise shows the potential of POMIC to use multiple patterns for selecting among multiple potential submodels and eventually select the most parsimonious and well fitting model version. We conclude that POMIC, although being a heuristically derived approach, can greatly improve the POM framework.     

A Bayesian approach to multi-source forest area estimation

In efforts such as land use change monitoring, carbon budgeting, and forecasting ecological conditions and timber supply, there is increasing demand for regional and national data layers depicting forest cover. These data layers must permit small area estimates of forest area and, most importantly, provide associated error estimates. This paper presents a model-based approach for coupling mid-resolution satellite imagery with plot-based forest inventory data to produce estimates of probability of forest and associated error at the pixel-level. The proposed Bayesian hierarchical model provides access to each pixel’s posterior predictive distribution allowing for a highly flexible analysis of pixel and multi-pixel areas of interest. The paper presents a trial using multiple dates of Landsat imagery and USDA Forest Service Forest Inventory and Analysis plot data. The results describe the spatial dependence structure within the trial site, provide pixel and multi-pixel summaries of probability of forest land use, and explore discretization schemes of the posterior predictive distributions to forest and non-forest classes. Model prediction results of a holdout set analysis suggest the proposed model provides high classification accuracy, 88%, for the trial site.

基于Monte Carlo法的城市雨水径流污染负荷模拟

为了评价流域雨水径流对受纳水体的污染程度和为城市雨水管理、受纳水体水质控制及规划提供依据,把以物理过程为基础的模拟特征融入到概率模型中,建立了城市雨水径流污染负荷的概率模型,并用Monte Carlo随机抽样法,对长期的径流污染负荷进行了模拟试验,得出流域径流污染负荷的分布特性及其特征值。以某市一小区的TP负荷为例进行了验证,模拟结果与实测结果较为一致。     

Does the African social spider Stegodyphus dumicola control the sex of individual offspring?

By scoring the chromosome number of developing embryos, we show that the sex ratio bias of the African social spider Stegodyphus dumicola Pocock is the result of an overproduction of female embryos. Only 17% of 585 embryos sexed from 14 egg sacs were male, a significant departure from a 1:1 sex ratio. We also explored the possibility of direct control of the sex of individual offspring in this species by examining the variance in the number of males per sac and the spatial distribution of male and female embryos within the sacs. We postulated that a variance in the number of males per sac lower than binomial (i.e., underdispersed or precise sex ratios) or a non-random distribution of male embryos within the sacs would suggest direct control of the sex of individual offspring. We found that the variance in the number of males per sac was indistinguishable from binomial and significantly larger than expected under exact ratios. Likewise, the spatial distribution of male embryos within three sacs examined was no more clustered than expected by chance. The sex ratio biasing mechanism in this species, therefore, apparently only allows control of the mean sex ratio but not of its variance. We present randomization and Monte Carlo methods that can be applied to test for departures from a random spatial arrangement of male and female embryos in an egg mass and for departures from binomial or exact ratios when not all members of a clutch have been sexed. Received: 21 October 1998 / Received in revised form: 23 March 1999 / Accepted: 26 April 1999     

Three long-range transport models compared to the ETEX experiment: a performance study

For operational or research purposes (dispersion computations of radioactive effluents during nuclear emergency situations, simulations of chemical pollution in the vicinity of thermal power plants), different models of passive dispersion in the atmosphere have been developed at the Environment Department of EDF’s R and D Division. This report presents the comparison of the performances of three such models: DIFTRA (lagrangian puff model, with operational goal), DIFEUL (three dimensional eulerian) and DIFPAR (Monte Carlo particle model) for the simulation of the first ETEX release, an international tracer campaign during which a passive tracer cloud has been followed over Europe. The results obtained in this study give model vs. experience differences of the same order as the model vs. experience differences observed during an international model comparison experiment using data of the Chernobyl release, the ATMES exercise. In addition to the standard statistical scores used in the evaluation of the performances of the transport models two asymmetric scores (in contradistinction with the Figure of Merit in Space) are proposed: “efficiency” and “power”. Their aim is to separate the two manners in which a model may be wrong: by predicting presence of pollutant while none is measured or conversely predicting absence when pollutant is actually detected.     

基于ANSYS-PDS的含缺陷压力容器可靠性分析

为更准确地分析含缺陷压力容器的可靠性,提出采用ANSYS软件中的PDS模块建模求解,从本质上克服API 581压力容器可靠性分析的局限性,采用将压力容器的内径、原始壁厚、缺陷深度、压力容器内压设置为服从正态分布的随机变量的方法,基于响应面法结合蒙特卡洛法抽样分析,通过编写APDL代码,探讨含缺陷压力容器的可靠性.结果表...